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In this lecture we discuss various aspects of QCD at nonzero chemical potential, 
^— H including its phase diagram and the Dirax; spectrum, and summarize what 

^ chiral random matrix theory has contributed to this subject. To illustrate the 

importance of the phase of the fermion determinant, we particularly highlight 
the differences between QCD and phase quenched QCD. 
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00 1. Introduction 



QCD at nonzero baryon chemical potential has turned out to be particularly 
challenging, and as of today, the phase diagram of QCD in the chemical 
potential temperature plane is fax from being understood. There is not 
even agreement of its gross features such as the existence of a critical end- 
point [1 3]. The reason for these difficulties is that the Dirac operator at 
nonzero chemical potential is nonhermitean resulting in a complex deter- 
minant so that the partition function cannot be evaluated by Monte-Carlo 
simulations. In this lecture wc compare the QCD partition function to the 
QCD like partition function that only differs by the absence of the phase of 
the fermion determinant. This theory, known as phase quenched QCD, can 
be studied by means of lattice QCD simulations and has a phase diagram 
that is certainly completely different from QCD. This makes it clear that 
the phase of the fermion determinant is absolutely essential for the physics 
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of QCD at nonzero chemical potential. 

In the second half of this lecture we will discuss various aspects of the 
phase of the fermion determinant and analyze its connections with the spec- 
trum of the Dirac operator. We will do this in the context of chiral random 
matrix theory [4,5] at nonzero chemical potential [6-8]. In the microscopic 
domain [5], this theory is equivalent to QCD and, in spite of the fact that 
it can be solved analytically, it shares the sign problem and the physical 
effects of the phase of the fermion determinant with QCD. 

2. Phase Diagram of QCD 

The phase diagram of QCD in the temperature-chemical potential plane 
can only be obtained by means of nonperturbative methods. The method 
that has had the most impact is lattice QCD, but because it requires a 
real action, only the temperature axis has been studied extensively. The 
consensus is that, in the temperature range of 150-200 MeV, QCD un- 
dergoes a crossover transition from a chirally broken phase to a chirally 
symmetric phase. The picture at nonzero chemical potential is much less 
clear. Although significant progress has been made for small chemical po- 
tentials, questions about the structure of the phase diagram along the chem- 
ical potential axis could not be answered based on first principles (except 
for asymptotically large chemical potentials [9]). That is why most of our 
knowledge of QCD at nonzero baryon density is based on models and gen- 
eral properties of the phases of QCD. 



Fig. 1. The phase diagram of QCD in the temperature baryon chemical potential plane 
(left) and in the temperature isospin chemical potential plane (right). The data points 
in the right figure are from lattice QCD simulations [10]. 

A tentative phase diagram of QCD in the chemical potential tempera- 
ture plane is shown in Fig. [l] (left). Among others we conjecture a chirally 
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symmetric crystalline phase made out of i? = ^ objects [11]. The best evi- 
dence for the existence of such phase is based on the Skyrme model where 
this phase exists as a strongly bound crystalline state with a binding en- 
ergy of 135 MeV per nucleon (see Fig. |2] (left)). The chiral restoration phase 
transition to this phase is of second order for a face centered cubic crys- 
tal [12,13] but is of first order for a cubic crystal [14]. The restoration of 
chiral symmetry in the dense phase is illustrated by the vanishing of the 
average cr- field over a unit cell (see Fig. [2] (right)). When the temperature 
is increased, we expect that this crystal will melt at a temperature that is 
of the order of the binding energy. This phase has also been identified as 
the quarkyonic phase [15,16]. 




Fig. 2. Energy of a unit cell of a Skyrmion crystal (left) for a cubic crystal with cubic 
boundary conditions (obtained in [14]) and for a bcc and fee crystal (obtained in [12]) 
versus the size of the unit cell. In the right figure we show the average value of the cr-field 
for an fee crystal [12]. 



Our aim is to understand QCD at ^ 7^ from first principles, i.e. starting 
from the partition function 

([[detiD + mf + fi-fo)), (1) 
/ 

but the nonhermiticity of D + i^jq makes this a very challenging problem 
indeed. Let us first discuss the effect of the phase of the fermion determinant 
by studying the theory where this phase has been quenched. For Nf = 2 
this partition function, known as phase quenched QCD, is given by 

^IQCDi {\det{D + m + fijo)]^) = (det(D -I- to -I- /170) det(i:) + m - /X70)) 

(2) 
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Therefore, the chemical potential of phase quenched QCD can be inter- 
preted as an isospin chemical potential [17], and at low enough temper- 
atures, pions will Bose condense for /i > m-„/2 [18-20]. This partition 
function can be simulated by Monte-Carlo methods [10,21-24] (see Fig.[T|) 
and its phases are in agreement with theoretical expectations based on a 
mean field treatment of a chiral Lagrangian [18,25,26]. 

We conclude that for /i > the phase factor quark determinant 

completely changes the phase diagram. The would be pion condensate is 
nullified by this phase factor after averaging. To better understand the 
effect of nonhermiticity, we will analyze the QCD partition function in the 
microscopic domain where it is equivalent to a chiral random matrix theory. 
This theory is analytically solvable in the nonperturbative domain of QCD, 
but the sign problem has all complications of QCD. 



3. Random Matrix Model at /x ^ 

A chiral Random Matrix Theory (chRMT) is obtained by replacing the 
matrix elements of the Dirac operator by Gaussian random numbers [4,5] 

(3) 



iW + ^1 

with W a complex n x (n+v) matrix, and is a multiple of the identity [6] 
or an arbitrary random matrix [8]. The random matrix partition function 
in the sector of topological charge v is given by 

Z,{mf-^i) = {\ldet{D + mf)), (4) 
/ 

where the average is over the Gaussian distribution of W . This partition 
function has the global symmetries and transformation properties of QCD. 
At fixed 9 -angle it is given by 

Z(m/,0;M) =^AC(Ai)e^"'Z,(m/;Ai) (5) 

V 

with Mui^[i) s normalization constant that may depend on /i and v. 
In the microscopic domain of QCD, defined by 

"^'«-^, ■^A4cd»1 (6) 



the mass and the chemical potential dependence of the QCD partition func- 
tion is given by the random matrix partition function 
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The reason for this equivalence is that in this domain, because the Compton 
wave length of the Goldstone modes is much larger than the size of the box, 
both theories reduce to the same static chiral Lagrangian. The Goldstone 
mass may either refer to physical quarks masses or to complex parameters 
(ghost quark masses) to probe the Dirac spectrum. 

Since the kinetic term of the chiral Lagrangian does not contribute in 
a mean field analysis, chRMT correctly describes mean field results for the 
low-energy limit of QCD beyond the microscopic domain provided that 
m < Aqcd and fi < Aqcd- 



4. Lessons from Random Matrix Theory 

4.1. Lesson 1: Homogeneity of Dirac Eigenvalues 

At nonzero chemical potential, the Dirac operator is nonhermitean with 
eigenvalues that are scattered in the complex plane. On inspection of 
quenched lattice QCD Dirac spectra (see Fig. 3) two properties stand out: 
i) The Dirac spectrum has a sharp edge, and ii) the 
distribution of the eigenvalues is more or less homo- 
geneous. Both properties can be understood in terms 
of random matrix theory [6,7]. In essence they follow 
from the fact that eigenvalues of nonhermitean ran- 
dom matrices behave as repulsive electric charges in 
the plane. Because the chiral condensate is given by 

S(m)=lV^— , 
^ ^ V m + tXk 

it can be interpreted as the electric field at m of 
charges at . Therefore, a direct consequence of 
Fig. 3. Quenched Dirac the homogeneity of the Dirac spectrum is that the 

eigenvalues on a 4^ X 8 i i i • i i ; • i • i • ; i 

lattice [27] quenched chiral condensate increases linearly with 

m when m is inside the domain of the eigenvalues. 



4.2. Lesson 2: Quenched Limit 

The quenched limit is the limit where the fermion determinant is ignored 
in generating the statistical ensemble. For zero chemical potential this has 
been a reasonable approximation, but at nonzero chemical potential, the 
limit of no fermion determinant is not given by 

lim„^o((det(D + ^70 + to))"), 

but rather by 



September 25, 2008 21:12 WSPC - Proceedings Trim Size: 9in x 6in prmnOS-jv 



6 



lim„^o(| det(D + + m)|"), 
i.e. the quenched hmit of phase quenched QCD [6]. Although similar ideas 
were proposed before [28,29], this result was first demonstrated convincingly 
for random matrix theory where the averages can be evaluated analytically. 

Because phase quenched QCD is QCD at nonzero isospin chemical po- 
tential, a phase transition to a pion condensate occurs at ^ = m-j^jl for 
low temperatures. This phase is absent in full QCD where, at low temper- 
atures a phase transition to a phase with nonzero baryon density occurs at 
/i = TOAr/3. We conclude that quenching completely fails for /i 7^ 0. 



4.3. Lesson 3: Width of the Dirac Spectrum 



The phase quenched partition function undergoes a phase transition at 
/i = TO7r/2 . In terms of Dirac eigenvalues the only critical point occurs 

when the quark mass hits the do- 



quark mass m 




main of eigenvalues (remember that 
it has a sharp edge) . Therefore these 
two points have to coincide [28] . This 
leads to the following condition for 
the half-width rric of the Dirac spec- 
trum: 



2mcS 
~F2- 



= V 



(7) 



Fig. 4. Scatter plot of Dirac eigenvalues 
for quenched or phase quenched QCD at 
nonzero chemical potential. 



This result has also been derived 
from the static limit of the chiral La- 
grangian for phase quenched QCD 
[19] which contains two low-energy 
constants, S and F^, as parameters. 
Since, at nonzero chemical potential 
the low-lying Dirac spectrum is characterized by two quantities: the density 
of eigenvalues and the width of the spectrum, it is possible to extract S and 
from these variables. More sophisticated methods to extract these con- 
stants from the low lying Dirac spectrum have been proposed and have been 
successfully applied to lattice QCD [30-38]. For sufficiently small chemical 
potential, when perturbation theory applies, the width of the Dirac spec- 
trum increases linearly with /i rather than quadratically. 
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4.4. Lesson 4' Infrared Dominance 

The question we wish to address in this section is whether the fluctuations 
of the fermion determinant given by 

det(D + m + /ijo) = '[li>^k + m) (8) 

k 

are dominated by the infrared part of the Dirac spectrum. At = 0, because 

the low-energy hmit of QCD is given by ehiral perturbation theory for 
sufficiently smaU quark masses, this is apparently the case. This is less 
clear at ^ when the eigenvalues are complex. In particular, because the 
eigenvalue density for large A increases as VX^ one would expect that 
also the ultraviolet part of the Dirac spectrum contributes to the phase of 
the determinant. Although a nonzero chemical potential does not introduce 
new infinities once the theory is regularized at /U = 0, one cannot exclude 
finite contributions from large Dirac eigenvalues. 

When the chemical potential and the quark masses axe in the micro- 
scopic domain of QCD, though, we can use random matrix theory to show 
that the average determinant is determined by the small Dirac eigenval- 
ues. Further evidence for the infrared dominance of this quantity comes 
from lattice QCD at imaginary chemical potential, where a small number 
of low- lying eigenvalues reproduce the random matrix result [39]. 

Because of agreement with lattice QCD we know that c;hiral perturba- 
tion theory can be applied to phase quenched QCD when both the chemical 
potential and tjItt are of order of F^^ or less. This shows that the magnitude 
of the quark determinant is infrared dominated in this domain. Since the 
full QCD partition function in this domain does not depend on the chem- 
ical potential at low temperatures, in agreement with chiral perturbation 
theory, we conclude that the /i-dependence of the the phase of the fermion 
determinant resides in the infrared part of the Dirac spectrum. 



4.5. Lesson 5: Failure of Banks- Casher at /x 7^ 

The Banks-Casher relation [40] states that 

E = lim lim ' , (9) 

m^OV^oo V 

where p{\) is the density of Dirac eigenvalues and V is the volume of space- 
time. Although originally intended for a Hermitian Dirac operator, this 

relation correctly gives a vanishing chiral condensate for phase quenched 
QCD at /Lt ^ 0. However, for full QCD at /z ^ , the chiral condensate has 
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a discontinuity when the quark mass crosses the imaginary axis, but it does 
so without the occurrence of an accumulation of eigenvalues. The alternative 
mechanism that is at work has been understood in detail in random matrix 
theory [41,42]. For a comprehensive review we refer to the talk by SplittorfF 
also in this volume [43]. Below we will illustrate this mechanism for QCD 
in one dimension which can also be viewed as a random matrix model. 



A 




E(m) 






m 


• 

m 


\ 






^sinh \i 





Dirac spectrum of Id QCD 

Fig. 5. Schematic plot of eigenvalues of the Dirac operator for lattice QCD in one dimen- 
sion (left). The yellow dots denote the position of the eigenvalues for a single gauge field 
configuration, whereas the green ellipse shows the support of the spectrum in the ther- 
modynamic limit. The chiral condensate for one flavor versus the quark mass is shown 
in the right figure. 

The partition function of lattice QCD in one dimension is given by 



Z= / dUdetD, (10) 

Jueu(Na) 

where the integral is over the Haar-measure of U{Nc). The Dirac operator 
is given by the N x Nc matrix (for N lattice points) 



D = 



ml 







ml 







ml 



(11) 



\-ef'U/2 ■■■ -e"^ ml J 
The chiral condensate for one flavor is defined by 

^Efe Ads:nfc(Afe + TO) 



S(to) 



(12) 



Notice that the determinant has a complex phase. For U{1) the chiral con- 
densate can be evaluated analytically [44] with the result that is shown in 
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the right figure of Fig. 5. The amazing phenomenon, also known as the "Sil- 
ver Blaze Problem" [45], is that the chiral condensate is continuous when m 
crosses the ellipse of eigenvalues, but shows a discontinuity at m = where 
there are no eigenvalues. This can happen because the chiral condensate 
is determined by exponentially large (in the number of lattice points) con- 
tributions which cancel to give a finite result for the chiral condensate for 
-> oo. 



4.6. Lesson 6: Quark Mass and Average Phase Factor 

The severity of the sign problem can be measured through the expectation 
value of the average phase factor. A physical interpretation is obtained by 
defining the phase with respect to the phase quenched partition function: 

, 2i0x ^ (dct^(£' + m-H/i7o)) ^ Zqqd 
'^"^ {\det{D + m + ^l^oW)- Z^CiCul' 

At nonzero temperature the free energies of QCD and |QCD| are different 
so that the average phase factor vanishes in the thermodynamic limit. 

At zero temperature, the free energy of QCD and phase quenched QCD 
are the same for /i < 771^/2. Therefore, in the thermodynamic limit, the 
average phase factor is one for fi < At finite volume, a nontrival 

/^-dependence is obtained from 
the one-loop corrections which 
give rise to an 0(1) contribu- 
tion. 

In the microscopic limit the av- 
erage phase is given by the 
random matrix result that can 
be evaluated analytically (see 
Fig. 6 (right)). In the left part of 
the same figure we show a scat- 
ter plot of the Dirac eigenvalues 
of 40 200 X 200 random matrices. 
We conclude that, in the micro- 
scopic domain of QCD, the aver- 
age phase factor vanishes if the 
quark mass is inside the support 
of the Dirac spectrum. 
The average phase factor has been evaluated to one-loop order in chiral 
perturbation theory [46] , and the results are in good agreement with lattice 







i ' 
I" 


O.B - 


-I t - 




f 


-1 






/j,= 0-3 




T = 0,0 





Fig, 6, Scatter plot of Dirac eigenvalues for chi- 
ral random matrix theory at nonzero chemical 
potential (left) and the average phase factor as 
a function of the chemical potential (right) for 
m = 0.1, In the left figure the mass is indicated 
by a black dot. 
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QCD [47] . The mean field limit of this result agrees with random matrix 
theory in the limit where the microscopic variables are large. Not surpris- 
ingly, also in the domain of chiral perturbation theory, the sign problem is 
most severe when the quark mass is inside the domain of eigenvalues. 



4.7. Lesson 7: Distribution of Small Dirac Eigenvalues 

One of the greatest successes of chiral random matrix theory has been the 
exact description of the distribution of the low- lying eigenvalues of the Dirac 
operator [4,5,48] both at zero (see review [49]) and at nonzero chemical 
potential (see review [50]). In this section we show recent results for the 
Dirac spectrum at nonzero chemical potential and topology. Such lattice 
calculations have become possible because of the introduction of an overlap 
Dirac operator at nonzero chemical potential [51]. In Fig. 7, the radial 
distribution of the lowest Dirac eigenvalue is shown for different values of 
the topological charge in comparison with analytical results derived in [52] . 




Fig. 7. Radial distribution of smallest Dirac eigenvalue for fi in different topo- 
logical charge sectors using a Dirac operator that satisfies the Ginsparg- Wilson relation 
(histograms) [52] compared to the prediction of chiral random matrix theory (smooth 
curve) . 

4.8. Lesson 8: Equality Two Condensates 

The chiral condensate can be calculated in two ways 

- lim lim - /y , ^ YiiT, -ly -^—\ , 

m^o y^oo V \ ^ Afc + m / y^oo m^o V \ ^ \k + m / 

(13) 
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where S*^^^ is nonzero because of spontaneous symmetry breaking, and its 
value does not depend on the total topological charge. The second chiral 
condensate can be expressed as 

S(^)= lim l <nA..^oA.-).=i ^ (^4) 

The reason is that in the sector of topological charge v, there are v exact 
zero modes. The equality E^^-* = E^^-* requires a subtle reshuffling of the 
eigenvalues: For v = 1 the eigenvalues are shifted by approximately half 
a level spacing w.r.t. v = Q va order to satisfy the chiral Ward identity. 



Therefore, A^~^ « \j '^k~'^ '^'k+i- correctly normalize the ratio in (14 1 
we evaluate it for a finite Dirac operator with an iV x {N + v) nonzero 
off-diagonal block. Such Dirac matrix has v exact zero eigenvalues (and 
perhaps additional paired zero modes which, being a set of measure zero, 
can be safely ignored). As it stands, the ratio ( 14 1 is dimensionally incorrect. 
A dimensionally correct ratio is obtained by replacing the largest squared 
eigenvalue pair ioi v = Q by its square root, [AJ^p^]^ AJ^". For the chiral 
condensate we then obtain the approximate expression 

E(^) = (15) 

Using the Banks-Casher relation, the smallest nonzero eigenvalue is approx- 
imately given by 

A„,i„ « (16) 

resulting in E*^^^ « ^E^^' This calculation can be made rigorous (see [53]) 
with the result that both condensates become equal. Since we compare the 
ratio of two partition functions, it is essential that they are normalized cor- 
rectly. This has been studied in a random matrix framework [53] confirming 
the above results. 

Also within a random matrix framework, it turns out that for /i 7^ the 
two condensate become only equal after the partition functions have been 
correctly normalized [54]. 



4.9. Lesson 9: Test of Algorithms 

Because random matrix models are exactly solvable and show all essential 
features due to nonhermiticity, they are an ideal tool to test algorithms 
for lattice QCD at nonzero chemical potential. We mention two examples: 
i) The density of states method was analyzed and tested in a chiral ran- 
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dom matrix model at nonzero chemical potential [55] and has been applied 
successfully to lattice QCD [56,57]. ii) The radius of convergence has been 
determined [58] for algorithms that rely on Taylor expansion [47,59] or an- 
alytical continuation in /i [60,61]. 

5. Conclusions 

We have shown that the phase of the fermion determinant dramatically 
affects the physics of the QCD partition function. For example, the phase 
diagram of QCD and |QCD| are completely different. This is our main mo- 
tivation for studying the behavior of the average phase factor and try to 
understand its relation with the spectrum of the Dirac operator. In par- 
ticular, this has been done in the microscopic domain where QCD is given 
by a chiral random matrix theory that can be solved analytically. In this 
lecture we have discussed nine questions where random matrix theory has 
contributed significantly to our understanding of QCD at /i 7^ 0. Perhaps 
the most important contributions are the failure of the quenched approxi- 
mation and the explanation of the mechanism that results in a discontinuity 
of the chiral condensate when the quark mass crosses the imaginary axis. 
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